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Abstract. In view of the current discussion on the subject, an effort is made to 
show very accurately both analytically and numerically how the Drude dispersion 
model, assuming the relaxation is nonzero at zero temperature (which is the case 
when impurities are present), gives consistent results for the Casimir free energy at 
low temperatures. Specifically, we find that the free energy consists essentially of two 
terms, one leading term proportional to T 2 , and a next term proportional to T 5 / 2 . 
Both these terms give rise to zero Casimir entropy as T — > 0, thus in accordance with 
Nernst's theorem. 

PACS numbers: 05.30.-d, 42.50. Nn, 12.20. Ds, 65.40.Gr 

1. Introduction 



The thermodynamic consistency of the expression for the Casimir pressure at finite 
temperature T is of considerable current interest. The problem gets accentuated at low 
T, where according to Nernst's theorem S = —dF/dT — > when T — > 0. Here 5" is the 
entropy and F the free energy per unit surface area. We shall consider the standard 
Casimir configuration, namely two semi-infinite identical metallic media separated by 
a vacuum gap of width a. The media are assumed nonmagnetic with a frequency- 
dependent relative permittivity e(uj). The two surfaces lying at z = and z = a are 
taken to be perfectly planar and of infinite extent. A sketch of the setup is given in 
figure [H 

The present work is closely related to our recent paper p] in particular, and also 
to our earlier papers on the thermal Casimir effect [21 EJ HI El E] . 
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We start from the Lifshitz formula: 



1 00 /-co . , 

(3F = —Y j ' \n(l- Ae- 2qa ) + \n(l- Be~ 2qa ) qdq, (1) 

27T n J( m /c L 



m=0 



where 



A={^^ (TMmode) (2a) 
B=( S -^-\ (TE mode) (26) 



s + p 
2rrk 



mT, (3 = 1/kT (2 c) 



-1+p 2 , p=^. (2d) 

Here ( m are the Matsubara frequencies, s and p are the Lifshitz variables, and the prime 
on the summation sign means that the case m = is to be taken with half weight. 
The appropriate dispersion relation to use is the Drude relation 

*«> = 1 + (3) 

where u = i(, u p being the plasma frequency, and v the relaxation frequency. The 
plasma wavelength is X p = 2ttc/u p . Our motivation for adopting the form ([3]) is that it 
agrees well with permittivity measurements (performed at room temperature). In the 
case of gold, 

uj p = 9.03 eV, v = 34.5 meV, A p = 137.4 nm. (4) 

The Drude relation is good for (<2x 10 15 rad/s. For higher (, the relation gives too 
low values for the permittivity (cf. figure 1 in |6J). Actually, the numerical input data 
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we used were taken directly from tabulated data along the imaginary frequency axis, 
s(i() for C > (courtesy of Astrid Lambrecht). These data cover 7 decades: 

1.5 x 10 11 < C < 1.5 x 10 18 rad/s. (5) 

For £ < 1.5 x 10 11 values for e(i() are obtained from ([3]) by extrapolation, but by our 
numerical evaluations only the m = value fell within this region. 

As mentioned, the permittivity measurements are made at room temperature. For 
defmiteness, we shall in the following use the room-temperature value v = 34.5 meV 
already given in (Hj), although we expect that at very low temperatures the true value 
of v is actually lower - cf. the recent discussion on this point by Klimchitskaya and 
Mostepanenko [7j. This fact will change our results quantitatively, but not qualitatively. 
In particular, it will not change our main conclusion regarding the validity of the Nernst 
theorem when T — ► 0. 

Let us emphasize the main assumption underlying our calculations: We assume v 
to possesses a nonzero value, however small, at any fixed temperature including T = 0. 

The assumed constancy of v might be questioned, as the Bloch-Griineisen law 
predicts that v depends on T as (cf. Appendix D in [6]) 

u(T) oc T 5 , T -> 0. (6) 

Such a relationship is not followed in practice, however, since there are always impurities 
which give rise to nonzero resistivity and so nonzero relaxation frequency at zero 
temperature [5]. In practice, therefore, our assumption above is always satisfied. The 
important point is that the relationship 

C 2 [e(iC) - 1] - 0, C^O (7) 

is always satisfied. It implies that the zero-frequency TE mode does not contribute to 
the Casimir force. The first to emphasize this kind of behaviour were Bostrom and 
Sernelius [9j , and the issue was discussed in detail in |6] . There are several other papers 
arguing along similar lines. Thus Jancovici and Samaj pU] and Buenzli and Martin [TT] 
considered the classical plasma of free charges in the high-temperature limit, where only 
zero frequency contributes, and they found the linear dependence in T in the Casimir 
force to be reduced by a factor of 2 from the behaviour of an ideal metal (the IM model). 

To illustrate the magnitude of the Drude thermal correction to the Casimir pressure, 
we give in Table 1 some calculated values, in mPa. It should be noted that if T increases 
from 300 K to 350 K, we find that 

(i) if a = 0.2 jum, the Casimir pressure diminishes by 0.4%; 

(ii) if a = 2.0 /zm, the Casimir pressure diminishes by 3.7%. 

The optimum gap width in connection with Casimir thermal corrections thus seem 
to lie around a = 2 /xm. 

An argument that has been put forward against the Drude relation is that by 
omitting the zero frequency TE term one gets a term linear in T in the free energy. 
Such a term would lead to a finite entropy at T = and so come into conflict with 
Nernst 's theorem. There are several recent papers arguing along these lines, written 
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Table 1. The Casimir pressure (in mPa) between Au-Au plates versus gap width a, 
when T = {1, 300, 350} K. Data extracted from Ref. [3]. 



o//im 


T = 1 K 


T = 300 K 


T = 350 K 


0.2 


508.2 


497.8 


495.7 


0.5 


16.56 


15.49 


15.30 


1.0 


1.143 


0.9852 


0.9590 


2.0 


7.549 x 10" 2 


5.550 x 10~ 2 


5.344 x 10~ 2 


3.0 


1.520 x 10^2 


1.033 x 10~ 2 


1.049 x 10~ 2 


4.0 


4.858 x 10" 3 


3.481 x 10~ 3 


3.804 x 10" 3 



from somewhat different perspectives [7J [12l [T3J, [HJ ISl [TU [T7] , and it is argued there that 
for such lattices the Drude model violates the Nernst theorem. Perfect crystal lattices, 
without impurities at all, are considered in these papers. It is argued therein that for 
thermodynamical consistency, relaxation due to electron-phonon scattering present at 
finite temperature should be neglected, and the use of the plasma dispersion relation 

e(iC) = 1 + (8) 

or a generalized version thereof is presented. The relation ([8]) corresponds to setting 
v = in PJ), such that (JED does not satisfy the condition (J7j). The plasma relation 
leads to quite a small temperature dependence in the Casimir force (correction oc T 4 ) 
in contrast to the distinct and almost linear decay with the Drude relation. Actually, 
the Drude theory in the limit v — > preserves entropy S = at T = 0, but 5* changes 
more and more abruptly at T = the smaller v is. 

In the following we intend to show very accurately, both analytically and 
numerically, how the Drude relation with v ^ leads to results that are in full agreement 
with the Nernst theorem. 



2. Analytical approach 

We start from the Drude model assuming some constant value for z/, and consider in 
the following only Matsubara frequencies that are relatively small, ( (= ( m ) v. These 
frequencies are the crucial ones for the behavior in the T — > limit. It is always possible 
to consider these frequencies when z/, as mentioned above, is finite. Then, 

£ (iC) = u p ~E D = ^l (g) 
We consider only the TE mode, which is the mode of main interest. Replace q by x: 

* 2 =(F^ = g> C«, do) 

Then the TE mode coefficient ( l2fe|) becomes 

B = (VTTx* - x)\ (11) 
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and the TE part of the free energy can be written 



PF TE = Cj2>g(m), 

m=0 



where 



g(m) = m 



.x In 



1 — B exp ( \J D( x 



dx. 



Now invoke the Euler-Maclaurin formula: 



m=0 

One then finds that 



, 9(u)au-±g'(0) + ^g"'(0)- 



poo I 

g'{0) = / xln(l -B)dx = - T (21n2 - 1). 
Jo 4 

And thereby one gets 



AF 



TE 



c 



(2 In 2 - 1) 



1 u. 



^(fcT) 2 (21n2- 1) 



(12) 



(13) 



(14) 



(15) 



(16) 



48/T ' 48 c 2 ftz/ 

valid for T 0.01 K. This result was first given by Milton at the QFEXT03 workshop 
0- 

Including the leading correction (Euler-Maclaurin summation starting at m = 1 
instead of at zero), one gets [1] 



AF TE = - 



T2 m 



1 + 0.204^S + ... 



(17) 



W(0) 

For gold plates, with a = 1 fim 

AF TE = dT 2 [l - C 2 T^ 2 + ...], (18) 

with 

d = 5.81 x 10~ 13 (J/m 2 K 2 ), C 2 = 3.03 K" 1 / 2 . (19) 

In order to avoid negative values for T slightly larger than 0.1 K, it is convenient to 
introduce the Pade approximant form 

C{F 2 



th ~ 1 + CoT 1 / 2 ' 



(20) 



This is equivalent to ffTB"]) with respect to the first two terms. Results ffl8|) - (T2Q|) were 
first obtained in Ref. [I]. 



3. Numerical calculations 



In the numerical calculations we assume two gold plates, with a = 1 /im. All dispersive 
data needed are in the experimentally known region given by ([5]) above. As mentioned, 
the only place where there is a need to use Drude relation ([3]) explicitly, is when m = 0. 
Actually, it is immaterial whether we use the experimental Lambrecht data (JSj) or the 
Drude relation directly. Thus figure H] is calculated with the use of the Drude relation 
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for all frequencies, but it turns out that a practically identical figure (with some noise 
because of lower accuracy) is obtained if we use Lambrecht's data. 

At T = the free energy is calculated numerically as a double integral rather than 
a sum of integrals, using a two-dimensional version of Simpson's method with adaptive 
quadrature. 

As for the TM mode, it is known that for ideal or nonideal metals the temperature 
correction for this mode behaves as T 4 . Thus, it is a smaller correction than the T 2 and 
T 5 / 2 corrections associated with the TE mode. We repeat that the dependence of v on 
temperature is neglected, and that we employ the room-temperature values for v given 
in 05). 

The vanishing of the zero-frequency mode is connected with the behaviour of the 
coefficient B at vanishing frequency. To illuminate this point, we show in figure [2] both 
coefficients A and B as a function of imaginary frequency and transverse momentum k± 
for an interface between gold and vacuum. In part (c) of the figure, we see how B — > 
when ( — > for k± ^ 0, whereas A in figure [2(a) for the TM mode equals 1 for all k± 
when ( — >• 0. 

By direct numerical integration and lengthy summations independent of the 
analytical derivations made in the previous section, we obtain the free energy 
numerically. Figure [3] shows the free energy versus temperature up to 800 K, while 
the inset shows details of the parabolic shape close to T = 0. The figure shows the 
decrease of the magnitude of the free energy and thus also the related decrease of the 
Casimir force up to a certain temperature. The inset shows how the slope is horizontal 
at T = 0, as predicted. Thus the entropy at T = is zero, in accordance with Nernst's 
theorem. 



4. A more accurate test 



Now, there are always uncertainties connected with numerical calculations. It is possible 
to make a much more accurate and sensitive test of the behaviour near T = in the 
following way. Define the quantity R as the relative difference between the temperature- 
dependent theoretical free energy AF™, and the temperature-dependent numerical free 
energy AF I TE ■ 



imm' 



AF™ - AF TE 

tn num 



r> th num /tii \ 

R= Af,I E <21) 
Assume for AF™ the Pade approximant form (l20l) . and assume for AF^, the expansion 

AF n ™ = D 1 (T 2 - D 2 T^ 2 + D 3 T 3 + ...) (22) 

with calculated values for the coefficients Di,D 2 and D 3 . Then, 

R = + ^{D 2 - C 2 )T^ + ^(C 2 D 2 - D 3 )T + ... (23) 

If Ci = D x and C 2 = D 2 : 

R(T = 0) = 0, RocT, T -> 0. (24) 
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Figure 2. Squared reflection coefficients A and B of the metal interfaces for the TM 
and TE modes, as a function of £/c and the transverse momentum k±. a) A for the 
TM mode, b) B for the TE mode, c) B for k± and £ close to zero. 



Calculated values of R are plotted in figure HI We see that R, when extrapolated, 
approaches zero linearly with a finite slope. This demonstrates the accuracy of the T 2 
and T 5//2 terms in the free energy. 
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Plate material : Gold 




Temperature (K) 



Figure 3. Numerical evaluation of the free energy (p} between two gold half-spaces 
as a function of temperature. The inset gives details for low T. 




5. Alternative derivation by expansion of g(m) 

It may be of interest to mention that, as a variant of the analytic approach, the 
dependence of the free energy on T near T = can be found by means of complex 
integration. Start from the TE expression 

°° roc 

p F = Cj2' /. x\n(l- Be' ax )dx, (25) 



m=0 



>t/D 
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where 



(3hvc 

and expand the logarithm 



C = a = 2aV2TrCm, (26) 



! 

OO 



(3F TE = -Cj2™ x ]T B n e~ nax dx. (27) 

m=l n=l 

Now use the formula 



| rc+ioo 

— / ds(nax)~T(s), 4 > c> (21 
2m Jc-ioo 



and sum over m, 



E ^' S/2 = C(|-l)- (29) 

m=l VZ 7 

Here £ denotes the Riemann zeta function. Distorting the contour so as to encircle the 
poles of the T function at s — 0, —1, . . . then yields the same result (jTSJ) as above. 



6. Summary and further remarks 

The main point in our analysis has been to show both analytically and numerically that 
the Drude dispersion relation ([3]) does not run into conflict with basic thermodynamics, 
as long as v ^ at T = 0. As we have seen, quite an accurate analysis is needed for 
this purpose. If we instead had argued in a more crude way, simply setting the TE 
coefficient B m = for m = and keeping all the other coefficients A m and B m equal 
to 1 as in the modified ideal metal model (MIM), then we would have broken Nernst's 
theorem. This issue has been discussed at length in Refs. [6] and [IB] . 

Whether the Drude predictions for the Casimir force are correct or not is to 
be decided upon from experiments. A difficulty here is the inherent uncertainty of 
theoretical predictions due to the relatively large spread of published data for the 
dielectric permittivity for typical metals such as Au - cf., for instance, the recent 
discussions on this point by Pirozhenko et al. [19] and Munday and Capasso [20J. 
The experiment with the highest precision [T5J [16] apparently is in disagreement with 
the Drude model, or any model satisfying (GO). It has also been suggested that there are 
large thermal effects due to surface roughness [21J. We might note that the 1% precision 
in the dynamic measurement made by the Purdue group [T5l [16] is not matched by 
the 3% accuracy of the very recent dynamic experiment reported in Ref. [22]. Our 
main concern in the present paper, however, has been to discuss the consistency of this 
theory. We wish to point out that it would be quite strange if the Drude relation, 
proved to be representing permittivity measurements with great accuracy, should turn 
out to be inapplicable to explain Casimir force measurements. Let us also mention here 
that an interesting discussion about the thermal Casimir effect and the Johnson noise 
has recently been given by Bimonte [23], as a possible theoretical explanation for the 
discrepancy with experiment. 



Drude Model and Casimir Entropy 



10 



The basic assumption for our analysis ought to be re-emphasized. We assumed the 
relaxation frequency v to be a finite quantity, for any value of T. One might here raise 
the question: what happens if the metal is a perfect crystal, with no impurities at all? 
In such a case v{T = 0) = 0, and the formalism above becomes inapplicable. (In this we 
have an opinion different from the definitive claim of a violation of the Nernst theorem 
given in Refs. [T71 US], for example. See also Ref. [21] ■) On basis of the calculation 
above, we can thus make no firm statement about the validity of the Nernst theorem in 
this special case. 

We ought to mention, though, that on physical grounds there are conceptual 
difficulties in simply setting v = in the dispersion relation: 

(i) It would yield a contribution to the Casimir force from the zero frequency TE 
mode. This mode is however not a solution of Maxwell's equations and should therefore 
not occur. (A more detailed discussion can be found in Ref. [6], and in Sect. Ill in 
Ref. [25].) 

(ii) Introducing a zero TE mode for perfect crystals would imply that such a medium 
would behave differently from a real metal when taking the limit v — 0. This would 
create a discontinuity in behaviour that we find unphysical. 

There are additional physical effects that we have not taken into account above: 

(1) One such effect is spatial dispersion [26], implying that the wave vector k is 
present in the dispersion relation. Then e = e(u, k) would become finite for finite k. 
Only the special case £(0,0) would be infinite, and it would not appear natural that this 
"measure zero" case should yield a finite contribution to the Carimir foreefl 

(2) Another effect that could have been taken into account is the anomalous skin 
effect [28l [29] . This effect occurs when the mean free path in the metal becomes much 
larger than the field penetration depth near T = 0. Again, no contribution to the 
Casimir force is found from the zero TE mode, and the Nernst theorem is satisfied. 

Finally, we refer to the very recent microscopic theory of the Casimir force at large 
separations, i.e. the classical limit, using statistical mechanics [30J - cf. also |31j and 
further references therein. These authors make use of a joint functional representation 
of both matter and field, enabling them to integrate out the field degrees of freedom 
entirely. Important in our context is that they find the TE modes not to contribute in 
this regime, and that the Casimir surface pressure is 

C(3)kT , , 

This is precisely as predicted by the Drude model in the same limit. This conclusion is 
further supported by Svetovoy's recent demonstration [32] of the cancellation between 
TE evanescent wave (EW) and propagating wave (PW) contributions for large distances, 
yielding Eq. (|30|) . while at short distances the TE EW dominates for the force between 
two metal plates or between a metal plate and a dielectric plate, resulting in a linear 
temperature term in the force. 

I It could be mentioned here that the contrary view has been expressed in Ref. [S7J . 
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